Nonequilibrium spin glass dynamics from picoseconds to 0.1 seconds 
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We study numerically the nonequilibrium dynamics of the Ising Spin Glass, for a time that 
spans eleven orders of magnitude, thus approaching the experimentally relevant scale (i.e. seconds). 
We introduce novel analysis techniques that allow to compute the coherence length in a model- 
independent way. Besides, we present strong evidence for a replicon correlator and for overlap 
equivalence. The emerging picture is compatible with non-coarsening behavior. 

PACS numbers: 75.50.Lk, 75.40.Gb, 75.40.Mg 



Spin Glasses [ij (SG) exhibit remarkable features, in- 
cluding slow dynamics and a complex space of states: 
their understanding is a key problem in condensed- 
matter physics that enjoys a paradigmatic status because 
of its many applications to glassy behavior, optimization, 
biology, financial markets, social dynamics. 

Experiments on Spin Glasses [l|, y] focus on nonequi- 
librium dynamics. In the simplest experimental proto- 
col, isothermal aging hereafter, the SG is cooled as fast 
as possible to the working temperature below the critical 
one, T < Tc. It is let to equilibrate for a waiting time, 
tw. Its properties are probed at a later time, t + 1^,. The 
thermoremanent magnetization is found to be a function 
of t/tw, for 10^^ < t/U, < 10 and t^ in the range 50 s — 



10 s[3| (see, however, [4j). This lack of any characteristic 
time scale is named Full-Aging. Also the growing size 
of the coherent domains, the coherence-length, ^, can be 
measured [51, |6|. Two features emerge: (i) the lower T is, 
the slower the growth of ^(tw) and (ii) £ ~ 100 lattice 
spacings, even for T ^^ T^ and t^ ~ 10^ sl5|- 

The sluggish dynamics arises from a thermodynamic 
transition at Tj^ [sj, |9|. There is a sustained theoreti- 
cal controversy on the properties of the (unreachable in 
human times) equilibrium low temperature SG phase, 
which is nevertheless relevant to (basically nonequilib- 
rium) experiments [10|. The main scenarios are the 
droplets[llj, replica symmetry breaking (RSB)'l2l, and 



the intermediate Trivial-Non- Trivial (TNT) picturejl3|. 

Droplets expects two equilibrium states related by 
global spin reversal. The SG order parameter, the spin 
overlap q, takes only two values q = i^EA • In the RSB 
scenario an infinite number of pure states influence the 
dynamics [12, 1 14 llSl . so that all — gsA < <Z < <Zea are 
reachable. TNT[13| describes the SG phase similarly to 



an antiferromagnet with random boundary conditions: 
even if q behaves as for RSB systems, TNT agrees with 
droplets in the vanishing surface-to-volume ratio of the 
largest thermally activated spin domains (i.e. the link- 
overlap defined below takes a single value). 

Droplets isothermal aging[16| is that of a dis- 
guised ferromagnet [36j . A picture of isothermal aging 
emerges that applies to basically all coarsening systems: 
superuniversality[l6J. For T < T^ the dynamics con- 
sists in the growth of compact domains (inside which the 
spin overlap coherently takes one of its possible values 
q — dzqEA)- Time dependencies are entirely encoded in 
the growth law of these domains, C(^)- The antiferromag- 
net analogy suggests a similar TNT Aging behavior. 

Since in the RSB scenario q = equilibrium states do 
exist, the nonequilibrium dynamics starts, and remains 
forever, with a vanishing order parameter. The replicon, 
a critical mode analogous to magnons in Heisenberg fer- 
romagnets, is present for all T < T^[17]. Furthermore, q 
is not a privileged observable (overlap equivalence 14]): 
the link overlap displays equivalent Aging behavior. 

These theories need numerics to be quantitative [l^jllSl, 
M, [m, E 0, [13, S 111, but simulations are too short: 
one Monte Carlo Step (MCS) corresponds to 10"^^ s[l[. 
The experimental scale is at lO^'' MCS (~ 100 s), while 
typical nonequilibrium simulations reach ~ 10~^s. In 
fact, high-performance computers have been designed for 
SG simulations[23,[2i,[2i. 

Here we present the results of a large simulation cam- 
paign performed on the application-oriented Janus com- 
puter 29|. Janus allows us to simulate the SG instanta- 
neous quench protocol for 10^^ MCS (~0.1 s), enough to 
reach experimental times by mild extrapolations. Aging 
is investigated both as a function of time and temper- 




FIG. 1: (Color online) Fit parameters, A and a {C{t,tvj) — 
A{ty,){l + t/tw)"^^"^*™^ ) vs. i„ for temperatures below T^ 
{T = 0.6 line: fit (for i„ > 10^) to a(iw) = ao + "i logiw + 
a2log^tw, ao = 6.35795, qi = 0.18605, ^2 = -0.00351835, 
diagonal x^/dof = 66.26/63). Coherent oscillations are due to 
the strong correlations of a(tw) at neighboring times (neither 
statistical errors nor the fitting curve nor x'^/dof vary if one 
bins data in blocks of 5 consecutive tw). 




FIG. 2: (Color online) Left: SG coherence length ^i_2 vs. 
waiting time, for T<rc. Right: ^1,2 vs. h, (^1,2 oc ll^^^'"''). 
Also shown data for the 2-D site-diluted Ising model {L = 
4096, 25% dilution, average over 20 samples, r = 0.64rj''"'5, 
^1,2 and Ji rescaled by 2 for clarity). Full lines correspond 
to Ising a — 0, (coarsening) and to the SG, a(rc) = 0.616 301]. 
Inset: Kf',2(iw) - C?°2(tw)]/L vs. C?°2(tw)/L for r = 0.8 and 
L = 24,40 and 80 (Cm(U) from a fit ^i,2iU) = A{T)tl/''^'^'^ 
for L = 80 in the range 3 < ^1,2 < 10, see text). 



ature. We obtain model-independent determinations of 
the SG coherence length ^. Conclusive evidence is pre- 
sented for a critical correlator associated with the repli- 
con mode. We observe non trivial Aging in the link corre- 
lation (a nonequilibrium test of overlap equivalence [ij]). 
We conclude that, up to experimental scales, SG dynam- 
ics is not coarsening like. 

The D = 3 Edwards- Anderson Hamiltonian is 



n^- 



E 



Jx,yO'x <Jy , ((• . .) : nearest neighbors) . (1) 



The spins ax = ±l are placed on the nodes, x, of a cubic 
lattice of linear size L and periodic boundary conditions. 
The couplings Jx^y — ±1 are chosen randomly with 50% 
probability, and are quenched variables. For each choice 
of the couplings (one sample) , we simulate two indepen- 
dent systems, {cri; '} and {ax '}. We denote by (• • •) the 
average over the couplings. Model H]) undergoes a SG 
transition at Tc = 1.101(5)[33|. 

Our L — 80 systems evolve with a Heat-Bath 



dynamics [31|, which is in the Universality Class of the 
physical evolution. The fully disordered starting spin 
configurations are instantaneously placed at the work- 
ing temperature (96 samples at T = 0.8 « 0.73 Tc, 64 at 
T = 0.7 « 0.64 Tc and 96 at T = 0.6 « 0.54 Tc). We also 
perform shorter simulations (32 samples) at Tc, as well 
as i = 40 and T = 24 runs to check for Finite-Size effects. 
A crucial quantity in non equilibrium dynamics is the 
two-times correlation function (defined in terms of the 
&e\dcx{t,t^)=axit + t^)ax{t^))^,M^- 



C{t,t^)^L-3j2cx{t,t^)., 



(2) 



linearly related to the real part of the a.c. susceptibility 
at waiting time tw and frequency lo = ir/t. 



To check for Full- Aging [3] in a systematic way, we fit 
C{t,t„) as A(iw)(l + i/<w)"^/"^*"^ in the range t„ < 
t < lOtw (33, obtaining fair fits for all t„ > 10^ . To be 
consistent with the experimental claim of Full- Aging be- 
havior for 10^"* < t„ < 10"'^^ 3], a(tw) should be constant 
in this tw range. Although a(tw) keeps growing for our 
largest times (with the large errors injgZ] it seemed con- 
stant for tw > 10'*), its growth slows down. The behavior 
at iw = 10"'^^ seems beyond reasonable extrapolation. 

The coherence length is studied from the correlations 



(1) 



of the replica field qxitw) = o'x (iw)o'i; (t 



.(2)^ 



C4(r,iw) = T ^y^^qx{tw)qx+r{tv 



(3) 



For T < Tc, it is well described by[]J,l20| 

C4(r,tw)'-r-'^e-('-/«(*"»', a -0.5, 6-1.5. (4) 

The actual value of a is relevant. For coarsening dynam- 
ics a = 0, while in a RSB scenario a > and C4(r, i^) 
vanishes at long times for fixed r/^jty^). At Tc, the latest 
estimate is a = l + 7^ = 0.616(9) [30| . 

To study a independently of a particular Ansatz as ^ 
we consider the integrals 



hit^ 



drr''C4r,t^), 



(5) 



(e.g. the SG susceptibility is x^^(^w) = 47r/2(iw) )• As we 
assume L^^(t^) we safely reduce the upper limit to L/2. 
If a scaling form C4(r, tw) ~ ''""/(''/^(iw)) is adequate 
at large r, then //c(iw) oc [£,{tv,)]'''^^~"' ■ It follows that 
^k,k+iitw) = Ik+iitw)/Ikitw) is proportional to ^(iw) and 
Iiit^)<y^^tkli ■ We find ^(2)(i^)«o.8^i 2(iw), where ^(2) 
is the noisy second- moment estimate[9|. Furthermore, 
for Ci,2>3, we find Co,i(iw)«0.46a,2(tw), and S.^\t^) = 
l-066,2(tw), (C"' from a fit to (gD with a = 0.4). 



~:o.7 



(L=33*) 
'„ = 2''(i=20) 
'w = 2' (i=16) 
f„ = 2" (L=10 ) 




FIG. 3: (Color online) For appropriate iw and L, the nonequi- 
librium Ciink(i,iw) vs. C (t,iw) at T — 0.7 , coincides with 
equilibrium Qnnkl vs. q^ (full lines, data from [341]. see text 
for definitions). The length-time dictionary is L = 10 or 
iw = 2^^ L = 16 or iw = 2^^ and L = 20 or i„ = 2^^. 
The coherence lengths, ,t(2i2) ^2.75(3), C(2^^) =4.23(4) and 
^(2^3) = 5.40(7), are in the ratio 10:16:20. Hence, from^(2=*5). 
Fig. [21 we predict the equilibrium curve for L = 33 . 



Note that, when ^ -^ L, irrelevant distances r ;2> ^ 
largely increase statistical errors for Ik- Fortunately, the 
very same problem was encountered in the analysis of cor- 
related time series 32], and we may borrow the cure[38j. 

The largest t^ where L = 80 still represents L = oo 
physics follows from Finite Size Scaling|31|]: for a given 
numerical accuracy, one should have L > fc^i,2(tw) • To 
compute k, we compare ^f 2 for L = 24, 40 and 80 with ^^2 
estimated with the power law described below (Fig. ^— 
inset). It is clear that the safe range is L > 7^i.2(iw) at 
r = 0.8 (at Tc the safety bound is L > 6^i,2{tw))- 

Our results for ^1.2 are shown in Fig. [21 Note for 
T ~ 0.8 the Finite-Size change of regime at i^ — 10^ 
(6,2 -- 11). We find fair fits to C(iw) = A{T)tl^ '''•^^ : 
z{Tc) = 6.86(16), z(0.8) = 9.42(15), z(0.7) == 11.8(2) and 
z{0.6) = 14.1(3), in good agreement with previous numer- 
ical and experimental findings ziT) = z{T^) T^/T^ ^. 
We restricted the fitting range to 3 < ^ < 10, to avoid 
both Finite-Size and lattice discretization effects. Ex- 
trapolating to experimental times (tw^lO^** '^lOOs), we 
find ^=14.0(3), 21.2(6), 37.0(14) and 119(9) for r = 0.6, 
r = 0.7, r = 0.8 and T= l.l«Tc, respectively, which seem 
fairly sensible compared with experimental data[5!,[a]. 

In Fig. [21 we also explore the scaling of Ii as a func- 
tion of ^1^2 (-^1 oc ^^~°)- The nonequilibrium data for 
T — 1.1 scales with a — 0.585(12). The deviation from 
the equilibrium estimate a = 0.616(9) [SOJ is at the limit 
of statistical significance (if present, it would be due to 
scaling corrections). For T = 0.8,0.7 and 0.6, we find 
a = 0.442(11), 0.355(15) and 0.359(13) respectively (the 
residual T dependence is probably due to critical effects 
still felt at T = 0.8). Note that ground state computations 
for L < 14 yielded a(T = 0) w OA^. These numbers 
differ both from critical and coarsening dynamics (a = 0). 

We finally address the aging properties of Ciink(i,iw) 




Clink (t,tw) 



{x,y) 



c^{t,t^)cy{t,t^)l{iL^). (6) 



FIG. 4: (Color online) Bottom: Ciink(i,iw) vs. C^{t,t^) 
for T — 0.6 and some of our largest t„ (vertical line: gl^^ 
from[23|). We also show our extrapolation of the Cunk vs. 
C^ curve to iw = 10^* (~ 100s, see text). Top: Derivative 
of Clink with respect to C for T = 0.6 . The horizontal line 
corresponds to the slope of a linear fit of fw = 10^* extrapo- 
lations (the line width equals twice the error). Inset: As in 
bottom panel, for the ferromagnetic site-diluted D — 2 Ising 
model (same simulation of Fig. [2J. 



Experimentalists have yet to find a way to access Cunk, 
which is complementary to C{t,t^) (it does not vanish 
if the configurations at i -I- t^ and t„ differ by the spin 
inversion of a compact region of half the system size) . 

It is illuminating to eliminate t as independent variable 
in favor of C'^{t,tv,), Figs. [3] and [H Our expectation 
for a coarsening dynamics is that, for C^ < q^^ and 
large i^, Cunk will be C-independent (the relevant system 
excitations are the spin-reversal of compact droplets not 
affecting Cunk)- Conversely, in a RSB system new states 
are continuously found as time goes by, so we expect a 
non constant C^ dependence even if C < (7ea [3^ • 

General arguments tell us that the nonequilibrium 
Clink at finite times coincides with equilibrium correla- 
tion functions for systems of finite sue[lO|, Fig. [31 (Qunk 
is just C4{r — 1), while q is the spatial average of qx, 
Eq.iO). Therefore, see caption to Fig. [31 we predict the 
q^ dependency of the equilibrium, conditional expectation 
QiinkL for lattices as large as L = 33. 

As for the shape of the curve Cii„k = f{C'^^'tw)^ 
Fig. [31— bottom, the t^ dependency is residual. Within 
our time window, Ciink is not constant for C < qEA ■ For 
comparison (inset) we show the, qualitatively different, 
curves for a coarsening dynamics. Therefore, a major dif- 
ference between a coarsening and a SG dynamics is in the 
derivative dCunk/dC^ , for C^ < gl^. Fig. [H— top. We 
first smooth the curves by fitting Cunk = /(C^) to the 
lowest order polynomial that provides a fair fit (seventh 
order for t^ < 2^^, sixth for larger i^), whose derivative 
was taken afterwards (Jackknife's statistical errors). 

Furthermore, we have extrapolated both Ciink(^ = 
rtw,iw) and C(i = rt^,t^) to t^ « lO^^ (- 100s), for 



r = 8, 4, . . . , ^ [40| ■ The extrapolated points for tvj = 10^^ 
fall on a straight line whose slope is plotted in the upper 
panel (thick line). The derivative is non vanishing for 
C^ < q\p^ , for the experimental time scale. 

In summary, Janus [22| halves the (logarithmic) time- 
gap between simulations and nonequilibrium Spin Glass 
experiments. We analyzed the simplest temperature 
quench, finding numerical evidence for a non-coarsening 
dynamics, at least up to experimental times (see 
alsoi26j). Let us highlight: noneguj/i&rmm overlap equiv- 
alence (Figs. I3l4p : nonequilibrium scaling functions re- 
producing equilibrium conditional expectations in finite 
systems (Fig. |3|); and a nonequilibrium replicon expo- 
nent compatible with equilibrium computations [331] . The 
growth of the coherence length sensibly extrapolates to 
tw = 100 s (our analysis of dynamic heterogeneities [25l 
26| will appear elsewhere 35[). Exploring with Janus 
nonequilibrium dynamics up to the seconds scale will al- 
low the investigation of many intriguing experiments. 
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08533, FIS2007-60977, FPA2004-02602, TEG2007-64188; 
by GAM (Spain) and the Microsoft Prize 2007. 
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